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D Abstract 

We propose an integral formulation of the equations of motion of a large class of field 
theories which leads in a quite natural and direct way to the construction of conservation 

o 

connections, and its appropriate mathematical language is that of loop spaces. The 
equations of motion are written as the equality of an hyper-volume ordered integral to 

ON 

an hyper-surface ordered integral on the border of that hyper-volume. The approach 
applies to integrable field theories in (1 + 1) dimensions, Chern-Simons theories in (2 + 1) 
dimensions, and non-abelian gauge theories in (2 + 1) and (3 + 1) dimensions. The results 
presented in this paper are relevant for the understanding of global properties of those 
$_i theories. 



laws. The approach is based on generalized non-abelian Stokes theorems for p-form 
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1 Introduction 



Symmetries play a central role in the understanding of physical phenomena. The laws govern- 
ing the fundamental interactions in gauge theories and general relativity are strongly based 
on symmetry principles. On the other hand the developments of non-perturbative methods 
to study strongly coupled system rely on symmetries revealed by deep structures like the 
weak-strong coupling dualities in gauge theories. Even though the Noether symmetries of 
Lagrangians and equations of motion are very important in many aspects of a given theory, it 
is perhaps correct to say that the hidden symmetries are the ones that have proved to be most 
efficient in the development of exact methods for non-linear and non-perturbative phenomena. 
The best examples of that are low dimensional theories with applications in many areas of 
physics like condensed matter, integrable field theories and solitons. The hidden symmetries 
responsible for the solvability of those (1 + l)-dimensional theories appear in general as gauge 
symmetries of an auxiliary flat one-form connection A^. In fact, the connection is a functional 
of the physical fields, and the zero curvature condition for A^ is equivalent to the classical 
equations of motion of the theory. The crucial fact here is that the flatness condition imply 
that the path ordered integral of the connection between two given points is independent of 
the choice of the path joining them. That statement is a conservation law, and the conserved 
quantities are given by the eigenvalues of the operator obtained by the path ordered integral 
of A^ over the entire one dimensional space sub-manifold. The exact developments in soliton 
theories and in many non-linear phenomena in low dimensions, over the last decades, were a 
direct consequence of that very important and simple fact. 

If higher dimensional theories present similar structures or not is an open and interesting 
problem. There have been several approaches to tackle that question, and we want to discuss 
here one which is a quite straightforward generalization of the ideas described above. One 
should expect the conserved quantities of a d + 1 dimensional theory to be associated to 
integrals of quantities on the ci-dimensional space sub-manifold. However, that space can be 
seen as a path in a generalized loop space in the following way. Choose a reference point xr 
in the (d+ l)-dimensional space-time M, and defined the space of maps LM from the (d— 1)- 
dimensional sphere S' d ~ 1 into M, such that the north pole of S" 1 ' -1 is always mapped into xr, 
i.e. LM = {7 : S ,d-1 — > M | 7(0) = xr}. The images of those maps are (d — l)-dimensional 
closed hyper-surfaces £ in M based at xr, and each one of them corresponds to a point of 
LM. Given a rf-dimensional hyper-volume in M, one can scan it with a collection of those 
closed hyper-surfaces S. Such collection is a path in LM, and so the hyper-volume in M can 
be seen as a path in the loop space LM. The idea now is, for a given theory in M, to look 
for a one-form connection A in LM, such that the conditions for its curvature to vanish are 
equivalent to the classical equations of motion of that physical theory. The flatness condition 
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for A implies that its path ordered integral between two given points in LM (hyper-surfaces in 
M) is independent of the choice of path (hyper-volume in M) joining them. That would lead, 
in a similar way to (1 + l)-dimensional theories, to conserved quantities as the eigenvalues 
of the path ordered integral of A over the paths corresponding to the d- dimensional space 
sub-manifold. That is the approach put forward in [1] and implemented in several examples 
of field theories in (d + 1) dimensions. See [2] for a review of the interesting results obtained. 
Among the difficulties of the approach are those associated to the non-locality and to the 
reparameterization invariance of the physical quantities. Note that a given hyper-volume in 
M corresponds in fact to an infinite number of paths in LM, which is a consequence of the 
infinity of ways of scanning it with hyper-surfaces. So, the physical phenomena should not 
depend upon the change of scanning. Despite those difficulties it was possible to impose local 
conditions in M which lead to the vanishing of the curvature of the connection A in the loop 
space LM, and made the physical quantities reparameterization invariant [1, 2]. 

In this paper we want to use the very same ideas proposed in [1, 2] to construct conserved 
quantities for a large class of theories, which include integrable field theories in (1 + 1) dimen- 
sions, Chern-Simons theories with sources in (2 + 1) dimensions, and Yang-Mills theories in 
(2 + 1) and (3 + 1) dimensions. However, instead of looking for a connection in loop space 
which zero curvature condition is equivalent to the classical equations of motion, we propose 
an integral form of those equations, related to generalizations of the non-abelian Stokes the- 
orem, and which lead in a quite simple and direct way to the conservation laws. Consider a 
physical theory in a (d + l)-dimensional simply connected space-time M, and let Q be any 
(in some sense topologically trivial) <i-dimensional hyper-volume in M, and suppose that the 
dynamics of such theory can be described by integral equations of the form 

Pd _ ie f m r = p de f n j (L1) 

where dfl is the border of Q, and where Pd-i and Pd stand for hyper-surface and hyper-volume 
ordering integrations respectively. The quantities T and J are built out of d — 1 and d forms in 
M respectively, and which are functionals of the physical fields. The details of the construction 
will be given in the examples discussed in the next sections. However, we deal with local field 
theories and the equations (1.1) are a direct consequence of the local differential equations of 
motion of the theory and of some generalization of the non-abelian Stokes theorem. On the 
other hand, since (1.1) is valid on any hyper-volume Q, it turns out that (1.1) imply those 
local differential equations when Q is taken to be infinitesimally small. In order to define the 
ordered integrations in (1.1) we scan Q with (d — l)-dimensional closed hyper-surfaces based 
on a reference point xr on its border dVt. Therefore, the equations (1.1) are not really defined 
on each Q, but on the generalized loop space LQ = {7 : S^ 1 — > Q | 7(0) = xr}, i.e. the space 
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of mappings from the (d— l)-dimensional sphere S^ 1 to Q, such that its north pole is always 
mapped into Xr. Consequently, Q can be seen as a path in LQ, and the r.h.s. of (1.1) is 
defined on such a path, with the l.h.s. of (1.1) being evaluated on its end points. But there is 
an infinite number of paths in LQ corresponding to the same Q. When one changes the choice 
of path representing Q, both sides of (1.1) change. However, the non-abelian Stokes theorem 
leading to (1.1) guarantees that the changes are such that both sides of (1.1) remain equal. 
Therefore, (1.1) transforms "covariantly" under change of parameterization of Q. In addition, 
we show in the next sections that (1.1) transforms covariantly under gauge transformations 
associated to the differential forms leading to the quantities J 7 and J . 

An important consequence of the integral form of the equations of motion (1.1) is that if 
one considers a closed hyper-volume fl c , i.e. without border, then the l.h.s of (1.1) becomes 
trivial and one gets that 

P d e^c J = t (1.2) 

Note that Q c corresponds to a closed path in LQ C . Then let us choose an intermediate point on 
that path, i.e. a closed (d — l)-dimensional hyper-surface E, such that Q c = Qi + fl 2 , with Q\ 
being the part of Q c going from the infinitesimal hyper-surface E# around the reference point 
Xr to E, and Q 2 to the part going from E back to E#. Then, the ordered integration implies 
that Pd e Jn 2 P d e Jn i = t. The order may be reversed depending upon the definition of the 
ordered integration. By reverting the sense of integration along the path one gets the inverse 
operator, and so one can rewrite that relation as Pde n ^ J = Pde n ^ , with Vt^ 1 being the 
path Q 2 in reversed order. Since that is valid for any closed path passing through E# and E, 
one concludes that the operator P^e^^ is independent of the path Vt joining E# and E. That 
path independency is a conservation law, and by choosing appropriate boundary conditions 
as we explain in the next sections, one gets that the conserved charges are the eigenvalues 
of the operator obtained by the path ordered integral P^e^ 3 , with f2 corresponding to the 
whole space sub-manifold. In the examples we discuss such conserved charges are shown to 
be gauge invariant, and independent of the parameterization of the hyper- volumes as well as 
of the choice of the reference point xr. In the case where the space-time is of the form S x H, 
with IR being the time, and S being a space sub-manifold without border, i.e dS = 0, then 
one gets from (1.2) that P^e^s^ = fl. Therefore, such operator is not only constant in time 
but trivial. In many cases, that topological property of the space-time leads to quantization 
of charges. That is a very important consequence of our construction. 

Even though we have not introduced a one-form connection in the loop space LQ, that 
concept is hidden in the quantity J . In addition, since we use generalizations of the non- 
abelian Stokes theorem, the quantity J corresponds to some sort of curvature of the quantity 
J 7 , now seen as a connection on a lower loop space LdQ, made of the space of maps of the 
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(d — 2)-dimensional sphere S d ~ 2 to the border dQ of Q. Therefore, there must be some sort 
of Poincare lemma playing a role here, implying that the curvature of a connection which 
is already a curvature should vanish. So, in that sense J would play the role of the flat 
connection in the approach proposed in [1, 2]. Note however that the integral form of the 
equations of motion (1.1) does not require the introduction of a connection in loop space to 
obtain conserved quantities for the theories we consider in this paper. 

The paper is organized as follows: in section 2 we implement our construction for integrable 
field theories in 1 + 1 dimensions re-obtaining well known results in that research area using 
the integral form of the equations of motion (1.1). In section 3 we discuss the cases of the 
Chern-Simons theory in the presence of a source as well as the Yang-Mills theories, both in 
2 + 1 dimensions. An important result of this section is the quantization of the charges in the 
case where the two dimensional space sub-manifold has no border. In section 4 we discuss the 
interesting case of non-abelian gauge theories in 3 + 1 dimensions, in the presence of matter 
currents. An important result here is an integral formulation of the Yang-Mills equations, and 
the construction of gauge invariant conserved charges. In the appendices A and B we give the 
proofs of the non-abelian Stokes theorems used in our constructions. 



2 The case of curves: theories in 1 + 1 dimensions 

In a 1 + 1 dimensional space-time M we establish a dynamical equation relating a field g(x), 
element of a Lie group G, to another field C^(x), a 1-form taking values in the Lie algebra 
Q of G. The relation between those two fields is built as follows. Consider a path 7 in M, 
parametrized by a, and define a quantity W through the differential equation 

dW „ dx" TIr , . 

*F + C >*F W = (21) 



where x^, fi — 0,1, are the Cartesian coordinates in M of the points of 7. Integration of (2.1) 
can be formally written as W = P\e -'i M da ■ Wr, where Pi stands for the path-ordering, 
and Wr is an integration constant corresponding to the value of W at the initial point xr of 
the path 7. 

Given any smooth path 7 in M, with initial and final points denoted by xr and Xf 
respectively, we impose the following equation for the fields g(x) and C M (x) 

g(x f ) ■ g(xR)- 1 = P x e-Z-r^M^ (2.2) 

with g{xR) and g(xf) corresponding to the values of g(x) at the end points xr and Xf, and the 
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Figure 1: The two paths 71 = 7l • 70 and 72 = jt • 7-l, connecting the points x# and £/, 
which are used to construct the conserved charges as the eigenvalues of the operator (2.4). 
The horizontal paths are parallel to the space axis and the vertical ones to the time axis. 



r.h.s. of (2.2) is obtained by integrating (2.1) along the path 7, assuming that the integration 
constant is unit. Note that (2.2) has the form of (1.1) since the border of 7 corresponds to 
its end points, i.e. dj = and g{xf) ■ g(xpt)~ l stands for the integration of T on $7, 

which in this case would be a zero-form. 

The first important consequence of (2.2) is that the path ordered integral of C^(x) is 
independent of the path. Indeed, if 71 and 72 are two paths in M with the same end points 
xr and Xf, then (2.2) implies that Pi e~ ^ daC ^~d^~ = p 1 e ~ Jy 2 da ' c ^~fa r _ That fact together 
with some appropriate boundary conditions is a conservation law as we now explain. Consider 
the space-time as being M = K x R, and let x° = t and x 1 = x be the time and space 
coordinates respectively. We choose the coordinates of xr as being (t, x) = (0, —L) and of 
Xf as (t,x) = (t, L), with L being a length scale which will be taken to infinity at the end of 
calculations. We choose two paths joining xr and Xf as shown in Figure 1, i.e. the first path is 
7i = 7l - 7o an d the second 72 = Jt'l-L- Note that 74 and 70 are paths at constant time at t = t 
and t = respectively. On the other hand 7^ and 7_^ are paths at constant space at x = L 
and x = —L respectively. If one assumes that the time component of the one-form satisfies the 
boundary condition Co(t,—L) = Co(t,L), for all values of t, then the path ordered integrals 



of C M along 7L and 7_ L are the same, i.e. Pi e ^-l""^> ^ = p 1 e ~h L a ^^ x >^ = jj^y 
Therefore, the equality of the path ordered integrals of along 71 and 72 leads to the iso- 
spectral evolution equation 



daC^x)^- 



n -f daCJx)^- 



U(t) Pi e 



(2.3) 



Consequently the eigenvalues of the operator 

Q = P 1 e -X»^(*)^ = g( t ,L)g-\t,-L) 



(2.4) 



where in the last equality we have used (2.2), are constant in time. Similarly, one can express 
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those constants of motion as TrQ , for any integer N. When we take the limit L — > oo, one 
observes that the conserved charges are determined by the asymptotic values of the field g(x), 
which is a known result in soliton theory [3]. 

Let us now consider the case where the space-time is of the form M = S x H, where the 
space submanifold S has no border, like for instance the circle S = S 1 . It then follows from 
(2.2) that the path ordered integral of C M on the whole space S must be unity, i.e. 

since the initial and final points are the same and so <?(xr) = g{xf). That is true for any value 
of time, and consequently (2.5) can be interpreted as a conservation law, where the conserved 
charges are in fact trivial. However, depending upon the theory under consideration, one gets 
(topological) quantization conditions for some quantities. A simple example would be that of 
an abelian pure imaginary connection, = i J M , where (2.5) leads to JgdaJ^x)^ =2nn, 
with n integer. 

Note that (2.2) transforms covariantly under the gauge transformations 

C^hC^-d^hh- 1 g^hg (2.6) 

since (2.1) implies that P x e ~ L, d ° c ^ x )^r h(x f ) Pi e~ A daC ^ x ^ /i(x fi ) _1 , with x R and x f 
being the end points of 7. Therefore, the conserved charges given by the eigenvalues of (2.4) 
are invariant under those gauge transformations satisfying h(t, —L) = h(t,L). 

If 7 is taken to be a path infinitesimally short, so that its end points approach each other, 
then g at Xf can be written as an approximation of the value it has at xr by using a Taylor 
expansion: g(xf) = g{x^)+d^g{xR)8x^, with Sx 11 being the infinitesimal displacement between 
xj and xr. Therefore, the l.h.s. of (2.2), up to first order in 5x^, becomes g^x^g^XR)^ 1 ~ 
1 + d fl g(xn) gi^XR)' 1 5x^. In addition, the path-ordering effects in the integration of are 
of higher order in and therefore the r.h.s of (2.2), up to first order, becomes simply 
1 — C^Sx^. Since that is valid for any infinitesimal path located anywhere in the space-time 
M, we get that (2.2) implies the following differential equation for the fields g(x) and C M (x) 

C,(x) = -d,g(x)g-\x) (2.7) 

Therefore is of the form of a pure gauge field, and consequently its curvature vanishes, i.e. 

d,C u -d u C, + [C tl ,C u ] = (2.8) 

The relation (2.8) is the so-called Lax-Zakharov-Shabat equation [4] or the zero curvature 
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condition, which is the basic structure used in the development of exact methods in soliton 
theory and two dimensional integrable field theories. The equation (2.2) is therefore an integral 
formulation of the Lax-Zakharov-Shabat equation. One can in fact obtain one from the other. 
However, there are some subtleties in the integral formulation, since it works with the two 
fields and g(x) and a relation between them, namely (2.2). That approach leads in a quite 
natural way, as shown in (2.4), to the fact that the conserved charges come from boundary 
terms. Such result is known for a large class of soliton theories [3], but it is not so certain 
that it holds for integrable field theories not possessing solitons. It would be interesting to 
investigate that issue further. In addition, the integral formulation leads to the triviality of 
the charges, or its topological quantization, in the case where the space sub-manifold has no 
border. 



3 The case of surfaces: theories in 2 + 1 dimensions 

In the case of theories defined on a (2 + l)-dimensional space-time M the basic ingredient of 
our construction is the so-called non-abelian Stokes theorem for a one-form connection C^. 
Let £ be a two dimensional smooth surface on M, and let <9£ be its border, i.e. a closed curve 
on M. The theorem states that the path-ordered integral of C M around <9£ is equal to the 
surface ordered integral on S, of the curvature of C M , i.e. 

P ie -§ d ^c^ .W R = W R -P 2 e^ dTdaW ~ lG ^ w ^^ (3.1) 

where 

Gfiu — d^C v — d v C^ + [Cp , C u ] (3-2) 

A proof of (3.1) is given in the appendix A, but its meaning is the following. One chooses a 
reference point xr on the border of S and scan it with closed loops starting and ending at 
Xr. The loops are labelled by r such that r = corresponds to the infinitesimal loop around 
Xr, and t = 2tt corresponds to the border <9£. Each closed loop is parametrized by a such 
that a = and a = 2n corresponds to Xr. The l.h.s. of (3.1) is obtained by integrating the 
differential equation (2.1) along <9£, and Wr is the integration constant corresponding to the 
value of W at xr. The meaning of Pi is that such integration has to be path ordered. As 
shown in the appendix A, the r.h.s. of (3.1) is obtained by integrating on E the differential 
equation 

^ - VJ = (3.3) 
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with 

f 2n dr^ dr v 

J= / da W- X G, V W^-^- (3.4) 
J da dr 

and the meaning of P2 is that such integration has to be surface ordered according to the 
scanning of E with loops as explained above. Again Wr is the integration constant and 
corresponds to the value of V on the infinitesimal loop around Xr. That the two integration 
constants have to be the same can be understood by shrinking E to the reference point xr. 

We now show how to use the non-abelian Stokes theorem (3.1) to define an integral for- 
mulation of the Chern-Simons and Yang-Mills theories, both in the presence of sources, and 
on a space-time of 2 + 1 dimensions. We then show how to use (3.1) to construct conserved 
charges for those theories. 

3.1 Integral formulation of Chern-Simons theory with matter source 

Consider a theory on a (2 + 1) -dimensional space-time M for a vector field and a current 
J M , with /j, = 0, 1,2, and let its classical equations of motion be defined as follows. On any 
two dimensional smooth surface E on M, with border <9E, the fields must satisfy the integral 
equations 

where J^ u stands for the Hodge dual of the matter current i.e., J^ v = e^ up J p , and where e 
and k are coupling constants of the theory. The meaning of the path ordered (Pi) and surface 
ordered (P 2 ) integrals in (3.5) is the same as those in (3.1), i.e. the l.h.s. of (3.5) is obtained 
by integrating (2.1) with = ieA^, and its r.h.s. by integrating (3.3) with G^ v = — J^ v . 

Since (3.5) is valid on any E, then it has to hold true when E is taken to be an infinitesimal 
surface. It then follows that (3.5) implies local differential equations for the fields as we now 
explain. Indeed, take E to be a planar surface of rectangular shape on the plane defined by 
two axis of the Cartesian coordinates, let us say x^ and x", with \i and v fixed. The border <9E 
is then a rectangle of infinitesimal sides 8x^ and Sx u . We evaluate both sides of (3.5) by Taylor 
expanding the integrands around one given corner of the rectangle, and keeping things at the 
lowest non-trivial order. One can check that the l.h.s. of (3.5) becomes 1 + i e F^Sx^Sx" , 
with no sum in fi and z/, and where 

F ia , = d l A v -d v Ap + ie [Ap, A u ] (3.6) 

The r.h.s. of (3.5) in lowest order is given by 1+ — J pu 6x^5x u (no sum in \i and v). Therefore, 
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for an infinitesimal surface, (3.5) implies the local differential equations for 



F, 



[lis 




[IV 




1 



(3.7) 



which are the equations of motion of the Chern-Simons theory in the presence of an external 
current J M . 

On the other hand one observes that if one takes in (3.1) as = ieA^, and therefore 
Gpv = ieF^, and uses (3.7), then one obtains (3.5). In other words, (3.5) is a direct con- 
sequence of the non-abelian Stokes theorem (3.1) and the Chern-Simons equations of motion 
(3.7). Since (3.5) implies (3.7), we see that (3.5) is indeed an integral formulation of the 
Chern-Simons theory. 

Note that in obtaining (3.5) from the non-abelian Stokes theorem (3.1) we have dropped the 
integration constant Wr. That has to do with the covariance of (3.5) under gauge transforma- 
tions, as we now explain. The Chern-Simons equation (3.7) transforms covariantly under the 
gauge transformations ->■ g g~ x + \ d^g g' 1 , since -» g F^ g~ l , and g J^ u g~ l . 

From (2.1) we have that under a gauge transformation W — > gjW g" 1 , where gi and gf are 
the values of g at the initial and final points of the curve where W is defined. Consequently, 
on a closed curve one has that W c — > gnW c g^ 1 , where gn is the value of g at the reference 
point xr, where the curve starts and ends. In addition, J defined in (2.1), with G^ u re- 
placed by ^7 Jf+u, transforms as J — > gn Jg^, and so from (3.3) we have that V — > gnV g^ 1 . 
However, if W{ and V\ are solutions of (2.1) and (3.3) respectively, so are W% = W{k and 
V2 = hVi, with k and h constant group elements. Under a gauge transformation one would 
then have Wf — > gu W? g^ , and Vi — )■ gnVi g^ 1 , with i = 1,2. But since k and h are arbitrary 
group elements one should not expect them to depend upon A^, and so be insensitive to its 
gauge transformations. Therefore, one could as well conclude that W£ — > gn g^ 1 k, and 
V2 — > h gji V\ g^ 1 . The only way to establish a compatibility is to assume that k and h should 
belong to the center of the gauge group G, since gn can be any element of G. Since the inte- 
gration constants Wr in (3.1) have the same status in this discussion, as k and h, we have to 
take them to lie in the center of G to have the gauge covariance of the integral Chern-Simons 
equation (3.5). However, when that is done they drop out from (3.5) since they commute with 
the path and surface ordered integrals. So, when integrating (2.1) and (3.3) to construct the 
l.h.s. and r.h.s. respectively of (3.5) one should keep in mind that those operators can carry 
an integration constant lying in the center of G without destroying the gauge covariance of 
(3.5). That fact may be important in some applications. 

Note that both sides of (3.5) depend upon the choice of the reference point xr and also 
on the choice of the scanning of £ with loops. However, when one changes the scanning 
and the reference point, the non-abelian Stokes theorem (3.1) guarantees that both sides of 
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(3.5) change in a way that they remain equal to each other. In that sense one can say that 
(3.5) transforms "covariantly" under the change of scanning and reference point. In fact, even 
though we have defined the equation (3.5) on any surface E in M, it is formally defined on 
the loop space £E = {7 : S 1 — > E | north pole — > Xr G <9E}, consisting of maps from the 
circle S l into S, such that the north pole of S 1 is mapped into xr. The images of such maps 
are closed loops in E, starting and ending at xr. Since E is scanned by a collection of such 
loops, and the loops are points in £E, one can see E as a path in £E. Therefore, a change in 
the scanning of E corresponds to a change of path in £E representing the same physical E. 

Despite the fact that (3.5) is defined on loop space, it leads to very physical consequences, 
like conservation laws as we now explain. Consider the case where the surface E is a closed 
surface E c , i.e. with no border. Then, the l.h.s. of (3.5) is trivial and we are lead to 

P 2 e % &c d ^ d ° w~%,w^^ = j (3 8) 

Being a closed surface, E c corresponds to a closed path in the loop space £E C , starting and 
ending at the reference point xr. Consider now a point on that path corresponding to a loop 7 
in E c . It divide the path in two parts corresponding to two surfaces, i.e. we have E c = Ei + E 2 . 
From the ordering defined by (3.3) we have that 

P 2 e K JE 1 M dT P 2 e K JE 2 da dr — j (3.9) 

By reverting the sense of the integration along the path, one obtains the inverse operator when 
integrating (3.3). Then, Ei and E^ 1 are two surfaces corresponding to two paths starting and 
ending at the same points, namely the reference point xr and the loop 7, which is their 
common border. Therefore, (3.9) implies that the integration along two paths with the same 
end points gives the same operator. Since that is valid for any closed surface E c and any 
partition of it into two surfaces, we conclude that the quantity P 2 e^^ drd(jW 1 Jf> vW ib r it- 
is independent of the surface E. It depends only on the initial and final points of the path 
corresponding to E in loop space, i.e. the border <9E and the reference point xr on it. Note 
that such independency corresponds to the change of the reparameterization (scanning) of the 
surface, as well as to the change of the surface itself, but keeping the border and reference 
point fixed. Such surface independency leads to conservation laws as we now explain. 

First we consider the topology of space-time to be M = S x IR, with time being the real 
line H, and the space being the closed two dimensional surface with no boundary S. A simple 
case is when S is the two-sphere, i.e. S = S 2 . Then if we evaluate (3.8) in space (i.e., E c = S) 
we get that the quantity Qs = Pi e ^ § s dr da w lj '^ w ^ r ^ r = H ; is conserved in time and equal 
to unity. That is an interesting relation since it may imply that the net charge on the whole 
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space vanishes, or then that the charge must satisfy some quantization condition. 

For simplicity, let us consider the case of an abelian theory, with gauge group G = U(l), 
where the path and surface orderings are irrelevant. In such a case one has Qs = e^ g = 1, 
where q is the total charge in space, i.e. q = § s dr da J^v^r^r- Those equations establishes 
a quantization condition involving the total amount of charge in space and the Chern-Simons 
coupling constants, i.e. 

q 6 

— = 2 7rn with n integer (3.10) 

K 

Such result has two important consequences for Chern-Simons theory on the space-time M = 
5xR The equations of motion (3.7) imply that the time component J° of the current, 
namely the charge density p, is proportional to the space components of the field tensor, 
which is the pseudo-scalar magnetic field B, i.e. p = kB. So, the effect of the Chern-Simons 
equation of motion is to attach magnetic flux to the electric charge [5]. For a point particle 
we have p = e5(x — xi), with x,\ being the position vector of the particle. Therefore, the 
magnetic flux associated to such a particle is $ = -, and (3.10) implies it is quantized as 
$ = The second consequence of (3.8) and the topology M = S x H, is that the phase 
gained by a non-relativistic particle that moves around another, due to a Aharonov-Bohm 
type interaction, is no longer dependent on the Chern-Simons coupling constant k. For iV 
such particles the charge density reads p = e S(x — x a ) and the attached magnetic field 
B = - Y2a — Xg). After a double interchange of two particles, their phase exchange is 
given by A9 = ^— [5j. Due to the quantization condition (3.10) and using that q = Ne we 
get A6 = which is a rational number, and not any number as would be the case if (3.8), 
and so (3.10), was not used. 

Let us now consider the case where space-time is 1R 3 , and discuss how conserved charges 
can be constructed using (3.8). As we have seen, as a consequence of (3.8), the quantity 

7(E) =P 2 e^^rdaW-^ v w^^ (3n) 

is independent of the surface E, as long as its border <9E is kept fixed and also the reference 
point xr on it. In addition it is independent of the scanning (parameterization) of E with 
loops. We shall consider two surfaces, Ei and E 2 , with the same borders as follows. The first 
one, shown in figure 2, is made of two parts. The first part is a disk V { £ on the plane x l x 2 at 
time x° = 0, and of a radius which will be taken to be infinite. The second part is a cylinder 
x /, where J is a segment of the x°-axis going from x° = to x° = t, and is a circle 
of infinite radius parallel to the plane x 1 x 2 . We take the reference point xr to be on the 
border of the disk T>^?, as shown in 2, and scan Ei = X>£? U (S^ x /), with loops starting and 
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Surface S x = 2?&? U (S 1 ^ x /) Surface S 2 = (S* x l) U 

Figure 2: The surfaces Ei and E 2 , with the same border S^, , and reference point xr, used 
in the construction of conserved charges. 

ending at xr, labeled by r, such that for r G [0, it] we scan with r = corresponding 

to the infinitesimal loop around xr, and t = tt corresponding to S^, the border of T>^ . For 
t G [7r,2 7r] we scan x I with loops, which start at Xr, go up in the x° direction upon to 
x° = t' G /, go round S^, and come down to Xr again. By varying t' within I we scan the 
cylinder x I. Following the notation of (3.11), and the ordering defined by (3.3), one gets 

V(Ei) = V(2>£>) V(Sixl) (3.12) 

The second surface E 2 , as shown in figure 2, is also made of two parts. The first part is a 
cylinder Sq x /, with I being the same time interval as above, and S$ a circle parallel to the 
x 1 x 2 plane with infinitesimal radius. The reference point xr is on the border of the base of 
such cylinder at x° = 0. The second part is a disk T>j£ on the plane x 1 x 2 at time x° = t, and 
of a radius which will be taken to be infinite. The surface E 2 = (5*q x 7) U T>^ is scanned 
with loops starting and ending at xr, labelled by r, such that for r G [0, 7r], we scan Sj x J 
with loops, which start at xr, go up in the x° direction upon to x° = t' G I, go round Sq, and 
come down to Xr again. By varying t' within I we scan the cylinder Sq x I. For r G [tt,2 tt] 
we scan with loops which start at xr, go up to x° = t, go round a closed loop on T>^, 
and come down to xr again. By keeping the two legs going up and down fixed and varying 
the closed loops on T>®, we scan it entirely. Again, following the notation of (3.11), and the 
ordering defined by (3.3), one gets 

V (E 2 ) = V (So 1 x /) V (2>g) (3.13) 

Since Ei and E 2 have the same reference point and the same border, namely S^ at x° = t, 
we have from the surface independency of (3.11) that V(Ei) = V(E 2 ). We now impose the 
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following boundary condition on our system 



Jl2 = Jo ~ T ^ fOT r ^ °° ( 3 ' 14 ) 



r 



with 6 > 0, r 2 = (x 1 ) + (x 2 ) , and T (f) being an element of the Lie algebra of G, depending 
on the spatial direction defined by f = -. That condition implies that the quantity J = 
^ f Q 2w da W^J^W^r^, vanishes on loops at spatial infinity, and therefore from (3.3) 
one gets that V (S^ x I) — 1. In addition, since the circle Sq has vanishing radius we 
also get that V (Sq x I) = 1. Therefore, from the equality of (3.12) and (3.13) we get that 
V \T>®) = V (v!^) ■ Those two operators are calculated using the same reference point 

xr, which lies at the border of T>6$ . Consider now a reference point x® , which have the 
same space coordinates as xr, but at a time x° = t, i.e. it lies on the border of T>% }, 
just above xr (see figure 2). By changing the reference point the quantity J changes as 
J — > W~ x (x®, xr) J W(x®, Xr), where W(x® ,Xr) is obtained by integrating (2.1) on the 
path joining xr to x®. Therefore, if one now integrates (3.3) on T>® with this new reference 
point, one gets that V XR (V®) = W -1 (x^ , xr) V (*) (T>®) 1^(1^,3;;;). Therefore, one gets that 

X R 

Vjg (2>£) = W (x®, xr) V Xr W- 1 (x®, xr) (3.15) 

where the subindices indicate which reference point is being used in the integration of (3.3). 
Note that in this way V x (t) \ T>®) and V XR (v^) correspond to surface ordered integrals over 
the entire space at times x° = t and x° = respectively, and with reference points at spatial 
infinity (border of the infinite disks) and at the same times. Consequently (3.15) constitute 
an iso-spectral time evolution for the operator 

m = ^ e ^"^ 1J ^" = P ie -^^- (3.16) 



where in the last equality we have used the Chern-Simons integral equation (3.5), and where 
the spatial circle with infinite radius stands for the border of T>® . Therefore, its eigenval- 
ues, or equivalently Tr V x (t) (t>®) , are constant in time. Those are the conserved quan- 
tities for the Chern-Simons theory. Note that such conserved quantities are gauge invariant, 
since under a gauge transformation we have that V^(i) (v®) — > g (x® ) V^t) (v®) g ( x r) > 
where g [xr) is the element of the gauge group, performing the gauge transformation, eval- 



uated at the reference point x®. The conserved quantities are also independent of the way 
we scan T>®, since we have already shown above that the operators of the type (3.11) are 
scanning independent. In fact, it was that property that lead to the conservation laws. In 
addition, the conserved quantities are independent of the choice of the reference point on the 
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border of Doc '■ That has to with the fact that by changing the reference point, Vjt) [Vooj 
changes by conjugation by an element W obtained by integrating (2.1) along a path on the 
border joining the two reference points. Consequently, its eigenvalues are unchanged. 



3.2 Integral formulation of Yang-Mills theory with matter source 

We now consider another theory on a (2 + l)-dimensional space-time M, with the same field 
content, i.e. a vector field and a current J M , with /i = 0, 1, 2, and with its classical equations 
of motion being defined as follows. On any two dimensional smooth surface S on M, with 
border <9£, the fields must satisfy the integral equations 



Pie- 



fax da \ A » 



p e * e Je * d ° W-^F^-p J^+ie [ , F v ] ) W*£- *g- 



(3.17) 



where e is the coupling constant of the theory, (3 is a free parameter, J^ u is the Hodge dual of 

£fiu P J p , Ffi is the Hodge dual of the curvature of the connection, 
= d fJi A u — dpA^ + ie [A^, A v \. The meaning of the path ordered 
(Pi) and surface ordered (P2) integrals in (3.17) is the same as those in (3.1), i.e. the l.h.s. of 
(3.5) is obtained by integrating (2.1) with 



the matter current i.e., J Ml/ \ 
i e F = -f F up and F 



C, = ie [a^ + (3F^ 



and its r.h.s. by integrating (3.3) with 



(3.18) 



= i e ( - (3 J uv + ie (3 



F F 

/1 1 f 



(3.19) 



In order to obtain the corresponding local equations of motion we consider the integral 
equation (3.17) on an infinitesimal surface £ of a rectangular shape, on the plane defined by 
two axis of the Cartesian coordinates, let us say x p and x u , with /x and v fixed. The border 
<9£ is then the rectangle of infinitesimal sides 5x^ and 5x u . Evaluating the r.h.s. of (3.17) in 
lowest order, and Taylor expanding the integrand around one given corner of the rectangle, 
we get 1 + ie (f^ — (3 J^ u + ie fi 2 F^ , F v \ 5x p 5x u , with no sum in \i and v. Analogously, 
evaluating the l.h.s. of (3.17) in lowest order, and Taylor expanding around the same corner, 
one gets 11 + {d^C v — dyC^ + [C^ , C u ]) 5x p 5x v ', again with no sum in \i and u, and with C M 
given by (3.18), and so 



ie 



F^ + 13 (DpF v -D v FA +ie$ 



F F 



(3.20) 
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Therefore, equating both sides of (3.17), in lowest order, one gets 



DpFv - D U F^ = -J^ (3.21) 

where D^* = <9 M * +ie [Ap, *}. Taking the Hodge dual one gets the Yang-Mills equations in 
(2 + 1) dimensions in the presence of mater currents 

D u F Ufl = J M (3.22) 

Note that if one takes the non-abelian Stokes theorem (3.1) with the connection C M given 
by (3.18), and so its curvature given by (3.20), one gets the integral equation (3.17) by 
using the Yang-Mills equations (3.21) . Therefore, (3.17) is a direct consequence of the non- 
abelian Stokes theorem (3.1) and the Yang-Mills equations. In this sense, (3.17) is an integral 
formulation of the Yang-Mills theory in (2+1) dimensions in the presence of matter currents. 

We have put therefore the Yang-Mills theory in the same footing as the Chern-Simons 
theory in the presence of matter currents, with its integral equation being given by (3.5). Con- 
sequently, most of the results we obtained for Chern-Simons are also valid for the Yang-Mills 
using the same techniques. For instance, the integral equation (3.17) transform covariantly 
under the gauge transformations A^ — > g A^g^ 1 + - d^g and J M — >• g J^g^ 1 . In addition, 
it transforms covariantly under re-parameterization of the surface £ with loops, and change 
of the reference point, as explained in section 3.1. But the most important result following 
from (3.17) is that if E c is a closed surface with no border, then its l.h.s. is trivial and so 

P 2 e ie &c dT daW-^F^-pJ^+ieplF^F^W^^ = j (3^3) 

That is the equivalent for Yang-Mills of the equation (3.8) for Chern-Simons, and it leads to 
conservation laws. In particular, for a space-time of the form M = S x IR, with S being the 
space sub-manifold, assumed closed with no border, like for instance the two-sphere S 2 , one 
gets that Q s = P 2 & dTda ^{f^-p j^+ie^F^F^w^^ = ^ ig constant in time and 

equal to unity. Again, it implies that the total charge in space vanishes, or then it may lead 
to quantization conditions. In the case of an abelian gauge group, for instance U(l), one gets 
that 

2 7T Tl 

$ — q = for n integer (3.24) 

e 

where $ = § s dr da F^^-^r is the total magnetic flux (or magnetic charge), and q = 
<f s dr da J^ v ^^r, is the total electric charge in space. 

Again following the reasoning used in the case of the Chern-Simons theory in section 3.1, 
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we get that (3.23) implies that the quantity 



V(E) = P 2 e ie ^ dTd(TW ' 1 ^-^ J ^ +iep2 [ F ^ F A) w ^^ (3.25) 

is invariant under smooth deformations of the surface S as long as its boundary and reference 
point xr are kept fixed. In addition, it is also invariant under the change of the scanning 
of X with loops based at Xr. Those facts can be used to construct conserved charges for 
the Yang-Mills theory. Let us consider space-time to be H 3 , and let us assume that the 
space components of the field tensor and time component of the currents satisfy the boundary 
conditions 

Ju = Jo~ ^T(f) ^-^T'ir) for r ^ oo (3.26) 

with 5 , 5' > 0, r 2 = (x 1 ) 2 + (x 2 ) 2 , T (f) and T' (f) being elements of the Lie algebra of the 
gauge group G, depending on the spatial direction at infinity defined by f = r -. Consider now 
a disk T>% , , of infinite radius on the plane x 1 x 2 , at a given time t, and let it be scanned with 
closed loops starting and ending at a reference point x^ on its border. Consider the following 
operator obtained by integrating (3.3) on Poo, with G^ v given by (3.19) 

(3.27) 

where in the last equality we have used the integral equation (3.17), and where is the 

border of Pod 1 , i.e. a circle at spatial infinity. Then following the arguments used in section 

3.1, leading to (3.16), one concludes that the eigenvalues of the operator (3.27) are constant 

r ( (0\l w 

in time. Equivalently, one can write those conserved charges as Tr 7 m ( Poo ) • Again 
following those same arguments one concludes that such conserved charges are gauge invariant, 
and independent of the scanning of PS with loops, and also on the choice of the reference 
point x^ on its border. 



3.3 Zero curvature representation for gauge theories 

We now comment on the connection of our integral formulation of Chern-Simons and Yang- 
Mills theories in (2 + 1) dimensions and the approach of [1, 2] using flat connections on loop 
spaces. 

Given the Lie algebra-valued 1-form C = C^dx^ and the 2-form B = ^B^dx^ f\dx v on M, 
we construct a 1-form connection A in the loop space LM = {7 : S 1 — > M \ north pole — > Xr}, 
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as [1, 2] 

I>2tT I n 



A= I daW- L B^W—5x u (3.28) 







where 5 stands for the exterior derivative on the space of all parametrized loops with base 
point x R , with 5 2 = and Sx fJ, (a) A5x"(a') = -5x p (a')/\5x^(a). The curvature J 7 = 5A+AAA 
is given by 

i r 27T dx x 

+ 5 da £ da' Ua - a') [B% (x (a')) - G™ (x (a')) , B% (x (a))] 



~ 0(a' - a) [BZ (x (a)) - Gf v (x (a)) , flJJ (x (a'))] j %^^W) A 

where is the curvature of C^, i.e. G^ v = d^Cy — d u C^ + [C M ,CV], and W is constructed 
out of by integration of (2.1). 

In the case of Chern- Simons theory we consider 

Cn = ic An and Bn U = —Jn V . 

K 



Lemma 2.1 of [2] claims that if 



Bfj.ii G^ 



then J 7 = 0. Then, for our choice of C M and B^ given above we see that 

T = Chern-Simons equation is satisfied. 

For Yang-Mills theory we take 

C fM = ie(A lx + /3F^ and B„ u = i e (f^ - (3 + ie /3 2 F„ , F v ) . 

The condition B^ v — G^ = leads to the Yang-Mills equation D u F UfJ- = J M , and therefore 
gives the zero curvature representation of this theory in loop space. 

Therefore, for the cases of Chern-Simons and Yang-Mills theories in (2 + 1) dimensions, 
the integral formulation approach and zero curvature on loop space lead to the same results, 
and also to the same conserved charges. 
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4 The case of volumes: theories in 3 + 1 dimensions 



For 3+1 dimensional theories we consider the generalization of the non-abelian Stokes theorem 
for a 2- form connection proved in appendix B. Given a three dimensional volume Q and its 
two dimensional border dQ, a closed surface, the theorem relates the surface-ordered integral 
of the connection W~ l B pv W along dVL with the volume-ordered integral of 



K = I drV {l daW~ 1 (D p B pu + D p B up + D u B pp )W^^^ r + 



, dx^ dx v dx p 
da dr dC, 



2- 



da I da' [BZW-FZWtBZi*)] ^(<x')^>) x 



o jo 




in Q: 



V R P 2 6^n drdaW^B^W^ = p ^J n dCIC y R ^ 



In the next few lines we summarize the construction of the above equation. A reference point 
xr is defined on the border of Q, and around it we construct an infinitesimal volume, whose 
border is the infinitesimal closed surface Sr. We then scan Q with closed surfaces based at 
the reference point xr. In doing so, we define the parameter (, such that ( = stands for 
the infinitesimal surface and ( = 27r, for the boundary dQ. During this variation the 
quantities V are calculated for each surface, in each step, through equation 

dV 

"^--VA = (4.2) 
ar 

with A = Jq W daW^B^W^r^rda. The surface E is scanned with loops parametrized by 
a G [0, 2tt], starting and ending at the reference point, and the parameter r labels these loops. 
The Wilson lines W are calculated on the loops using equation (A.l). The l.h.s of (4.1) is 
therefore obtained after integrating (4.2) over the surface dfl, with the ordering given by the 
way we scan it with loops. The quantity Vr is the integration constant (the value of V for 
the infinitesimal surface £r). 

Once the surface £ is closed, V can also be obtained from the equation (as we show in the 
appendix) 

dV 

--CV = (4.3) 
and this result is expressed in the r.h.s of the theorem above. 
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We now proceed to show how to formulate an integral version of Yang-Mills theory through 
this non-abelian Stokes theorem, and also, how it leads to conserved charges. 



4.1 Integral formulation of Yang-Mills theory with source 

We consider the Yang-Mills theory in (3 + 1) dimensions for a gauge group G and in the 
presence of matter currents J M . The classical equations of motion are given by 



D v F vp = 



D v F Ufl = J" 



(4.4) 



where F^ v = d p A u — d v A p + ie [A^, A v \, and F pu is the Hodge dual of the field tensor, i.e., 

Fp,u = 2 ^ pvp\ • 

One can obtain an integral equation for the Yang-Mills theory using the Stokes theorem 



(4.1) as follows (see [6] for more details). Take B 



flU 



ie 



aF„ v + P F, 



with a and being 

arbitrary constants, and using Yang-Mills differential equations (4.4) to replace D p B pv + 
DpB up + D U B PP by (—ie/3 J^vx) ■ With that, the quantity /C above is now given by /C = 
J^drV JV-\ with 



J 



x 



" T j I ■ ojxr-1 7 Ti/ dx v dx X 

da { ie(3W 1 J pvp W — + 



da dr d( 



+ e z 



da' 



riW 
up 



pv 



dx K dx» ( dx p {a') dx u {a) dx p {a') dx u (a) 



da' da 



dr 



d( 



d( 



dr 



a 



x 



aF u „ + (3F, 



and V is constructed by integrating (4.2) with A = ie f Q n daW 1 
and where we have used the notation X w = W~ x X W. Then we have [6]: 



pv 



Ty dx^ dx" 
da dr ' 



(4.5) 



which is a direct consequence of the generalized non-abelian Stokes theorem (4.1) and the 
Yang-Mills equations (4.4). On the other hand, the integral equation (4.5) implies the differ- 
ential equations (4.4), as we now explain. 

Equation (4.5) is defined for an arbitrary volume Q, and in particular, for an infinitesimal 
one. Take Q as an infinitesimal cube of sides 8y v and 5z x , with the indices fixed (see 
Figure 3). We choose the reference point xr to be at one of the vertices, and when opposite 
surfaces are scanned with loops based on xr one has to pay special attention to the fact that 
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XR gx 



Figure 3: The scanning of an infinitesimal cube. 



the signs of the velocities in r-direction changes. 

Now, considering only first order contributions to equation (4.5), the integrand in its l.h.s 
can be evaluated at any point on the cube's face, since the differences will be of higher order, 
and therefore, evaluating it on the face 5z x = one gets 3 — ie aF^ v + f3F fll , Sx^Sy" '. For 

the face at xr + 5z the contribution comes from ie ( W~ x otF uv + [5F UU W) 5x^5y u . 

V L J / (x R +5z*) 

For this term we have to expand both the Wilson line and the field strength so that their values 
at xr + 5z are approximated to their values in Xr. The equation for the Wilson line gives for 
an infinitesimal variation a — > a + 5a along the z direction, W/ R+gg x\ ~ \ — ieA\ (xr) 5x x and 

jii/\-> It \~8z) = 1 fj, r \< [■; i 5 'A' 

Sx ti Sy u Sz x . Doing the same for the other two pairs of 



Taylor expanding the field strength in this direction, F fMU (xR + 6z) = F^ u (xr) + d\F^ v (xR)5z x 



we end up with ieD\ 
faces one gets 



OtF^jj ~\- /SF^jy 



J X R 



P 2 e-/an^K^ + ^]^^ w 1 + ie(Dx[aF flu + (3F^\ + cyclic perm.) x Jx»5y u 5z x . 

For the r.h.s of (4.5), considering only the first order contributions, we notice that the com- 
mutator term is of higher order with respect to the first term involving only the current, for 
it has one more integration along the loop. Then, up to lowest order 

P 3e fn^vjv-i Kl + ie {3j^ x s x »6y u 6z x 

and clearly equating the previous result with this one, we get the set of Yang-Mills equations 
(4.4), as the coefficients of the parameters a and /3. 

Comparing the integral equation (4.5) with (4.1) one notices that the integration constants 
Vr are missing. One has to keep in mind that while (4.1) is a mathematical relation, (4.5) is 
a physical equation, and therefore gauge covariance is important, and in order to guarantee 
that those integration constants must lie in the center of the gauge group, as we now explain 
in detail. Consider the gauge transformation of the Yang-Mills field — > gA^g" 1 + ^d^gg -1 , 



3 The minus sign is due to the choice of the direction of scanning. 
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which implies the following transformations of the field strength F^ u — > gF^g^ 1 , and the 
matter current transforms in the same way, J M — > gJ^g^ 1 . From (A.l), the Wilson line 
changes as W — > g$Wg~ l where g± and gj stand for the value the gauge group element takes 
at the points X{ and Xf respectively; for the reference point we denote as gn the value of 
the gauge group element there. Then it is direct to see (by replacing W , F^ v and by the 
respective gauge transformed quantities given before) that all the important quantities related 
to the integral Yang-Mills equation (4.5) (such as V, J ', K, and ^4) transform as V — > gnVg^ 1 , 
J -» 9RJ9R 1 , K> -» gak-gR 1 and A -> gaAgx 1 . 

The l.h.s of (4.5) comes from the l.h.s of the non-abelian Stokes theorem previously pre- 
sented, which, in turn, is the result of the integration of equation (4.2). We notice that if V 
is a solution of this equation, then V = kV is also a solution, where k is a constant element 
of the gauge group. Under a gauge transformation one gets V — > guVg^ = gakV 'g^ 1 . On 
the other hand, since k is an arbitrary constant, and so insensitive to transformations of the 
gauge field, one could have written it as V — > kg R Vg~^ : and the only way to guarantee the 
compatibility is to have g^k = kg R , i.e., to have k in the center Z(G) of the group. The same 
argument can be applied to the r.h.s of (4.5), which comes from the r.h.s of the non-abelian 
Stokes theorem, which, in turn, is the solution of equation (4.3). In that case, if V is a solu- 
tion, so is V = Vh, with h a constant element of G. Finally we conclude that (4.5) is gauge 
covariant only if the integration constants are in the center of the gauge group. Then they 
can be cancelled trivially, since they commute with the surface and volume integrals, and that 
is the reason why they do not appear in the integral Yang-Mills equation (4.5). 

An important issue concerns the fact that equation (4.5) is formulated in a way that it 
depends on the particular choice of the reference point xr and of the scanning of the volume 
with surfaces. Although it is at first sight unwanted, the Stokes theorem (4.1) guarantees that 
if one changes any of these things, each side of (4.5) will change in a way to remain equal 
to each other; in this sense, this equation transforms "covariantly" under reparametrization. 
This can be better understood once we realize that this equation is formulated in the loop 
space LQ = {7 : S 2 — > fi|north pole — > xr G dQ}, formed by maps from S 2 into Q, such that 
the north pole of S 2 is mapped into the reference point xr. The images of this map are closed 
surfaces in Q, starting and ending at xr. The volume Q is scanned by a family of these closed 
surfaces, which are points in LQ, thus, Q is a path in the loop space LQ. Then, a change in 
the scanning of Q corresponds to a change on the parametrization of the path in LQ, which 
does not change the path, nor the physical results from it. 

Equation (4.5) does not only describes the Yang-Mills theory in loop space, but also leads 
to a conservation laws as we now discuss. For a closed path in loop space, corresponding to a 
closed volume Q c in space-time, with no boundary, the l.h.s of equation (4.5) becomes trivial 
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and we get 

p^dSdrVJV-i =1 (46) 

Consider now a given point in that path, in space-time, the surface E, which divides it into 
two parts: Q c = Qi + Q 2 - Then, we can split equation (4.6) as 

Each term of this equation is obtained from integration of (4.3). It can be written as 
V(p,2) ' ^(^1) — 1- Now, reverting the order of integration for the second path O2, one 
gets equivalently V^ -1 ^^ 1 ) • V(Vli) = 1, which implies that V(£li) = V^fi^ 1 ). In other words, 
the integration of (4.3), for the volumes f2i and f^ 1 , with the same boundary £ and reference 
point xr, leads to the same operator. In fact, since equation (4.7) holds for any closed volume 
Q c and for any partition of it into two other volumes, what we just saw is that not only the 
integration of (4.3) along two volumes with same boundary gives the same operator but also 
that P 3 e^ d ^ dTV ^ v 1 is independent of the volume, for any volume fl This means that as long 
as the border is kept fixed, one can change the volume and there will be no consequences to 
V(Q). The fact that this operator depends only on the border dQ and on the reference point 
xr on it, leads to conservation laws as we now explain. 

Consider first the case in which space-time M has the topology S x IR, with S being a 
closed unbounded spatial submanifold and IR the time. Then, equation (4.6) can be evaluated 
for f2 c = S, giving that the quantity 

Qs^p^fs^yjy- 1 = i. (4.8) 

is conserved in time and is equal to the unit. This has two main implications: the first is 
that the net charge on the whole space vanishes, and the other is that it might lead to some 
quantization condition. In particular, for the case of the Maxwell theory, where the (abelian) 
gauge group is G = U(l), this equation for Q$ is satisfied if q = f s d(drdaJp U \^^^-, is 
such that 

2nn 



e/3 

with n integer. If for some reason can be fixed at the quantum level, then (4.8) express the 
quantization of electric charge [7]. 

Next we consider the space-time M to have the topology IR 3 x IR. Lets take two points 
in the loop space, which correspond in space-time to the following two surfaces (see figure 4): 
the infinitesimal 2-sphere Sq around xr, at time x° = and the 2-sphere which is the border 
of the entire space at time x° = t; and which we denote as . Now, these two points 
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51x7 



c2,(0) 



Figure 4: The schematic representation in loop space LQ 
of the "path" used in the construction of the conserved charge. 



in LQ are joined by two different paths (volumes in space-time) and each of these paths are 
composed by two parts as we now explain. The first part of the first path is the infinitesimal 
hyper-cylinder Sq x J, with / being the interval in 1R from to t and Sq the infinitesimal 

2 (t) 

sphere around xr. The second part of this path is the volume inside the sphere Soo , which 
we make by "blowing up" the infinitesimal sphere Sq when we reach the point x° = t; we 



denote this volume by tt®- This first path is then the composition fli = (Sq x I) U fiS. The 
second path has a first part, which in space-time corresponds to the whole spatial volume 



,(t) 



fiJx?, inside S^ yu ' at x° = 0, and its second part is the hyper-cylinder S'^ x /. This path is the 
composition f2 2 = ^Jx? x (<S^ x I). Basically we are dealing with different paths (volumes) 
which boundaries - Sq and S^^ - are common, and thus, following our previous result, the 
operator V(Q) = P 3 e^ d< ' dTVJV 1 , is independent of the volume, and should depend only on 
the boundaries of f2. Therefore, it follows that V(Q 2 ) = ^(^i), and this can be written as 



,2,(0) 



v(s<$ x j)v(n£)) = v(n£)v(s 2 Q x /). 



(4.9) 



In order to obtain V(Q) for each step it is necessary to evaluate fC = J Q 27T dr VJV~ l on the 
surfaces scanning each volume f2. We shall scan a hyper-cylinder S 2 x I with surfaces, based 
at xr, of the form given in figure (5.b), with t' denoting a time in the interval I. Each one 
of such surfaces are scanned with loops, labelled by r, in the following way. For < r < ^, 
we scan the infinitesimal cylinder as shown in figure (5. a), then for ^ < r < 4^ we scan the 
sphere S 2 as shown in figure (5.b), and finally for 4^ < r < 2n we go back to x R with loops 
as shown in figure (5.c). Then, we can split /C in three parts: K a + /C& + /C c , each of them 
corresponding to one of the three surfaces, defined by the r intervals. 

From the physical point of view it is very reasonable to take the current and the field 
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strength to satisfy the boundary conditions 



R 2+8 



for R — > oo, with 5 and 5' bigger than zero. With that, integration of J over vanishes 
and we get /Q, = for the path ^2- We notice that for the path fix, the infinitesimal sphere 
Sq does not contribute to /C&, and we have also /Q, = in this case. Then we conclude that K. 
calculated in both spheres and Sq gives the same result and therefore 



V(Sg x /) = V{Sl x /). 



It is now possible to contract the cylinders into a line, so that the loops scanning them in 
the intervals r G [0, ^] and r G [4^,27r] become exactly the same. However, since one set 
of loops is "going up", and the other is "going down" (the sign of the velocity ^ changes), 
the contributions K a and K c exactly cancel. Note that V inside /C, obtained from (4.2) does 
not change sign since it does not see the direction the loops are going to but only takes into 
account the profile of the loop for each r. 




2,(r 



(a) 



(b) 



(c) 



Figure 5: Surfaces of type (b) scan a hyper-cylinder S 2 x I. 

Finally what remains is exactly V calculated on fl^, whose border is the sphere 5» , 
where we need to evaluate /C. In order to scan the sphere, we establish the point , on 
its boundary. This is the point xr at time x° = t, so that the path starts at Xr, goes up 
to Xji, and from this point we scan the sphere with loops. Then, we go back to xr. So, 
we construct the Wilson line composed of the two parts: W = W(x, x^ )W(xr, Xr) the 
one corresponding to the "leg" that goes up from xr to X®, denoted by W(x^ ,Xr), and 
the one that corresponds to the path from x^ to the point x(a), on the sphere. Using this 
decomposition we can re- write every needed quantity in terms of the new reference point x^ . 
In particular, A in equation (4.2) is decomposed as A XR = W(xr, xr)~ 1 A (*) W(xr, xr). 
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This leads to V XR = W 1 (x$, Xr)V (*) W(x R , Xr); In the same way the quantities /C& and J 

X R 

transform as K b>XR = W~\x^, x R )IC bx (t)W(x R , x R ) and J XR = W~ 1 (x^, x R )J x <$ W(x R , x R ), 
so that V(Q) obtained from (4.3) becomes 

v XR (n£) = w-\ x { lx R )v XR (n£)w(x%,x R ). 

Plugging this into the equation (4.9) we finally get 

v* H (n®) = u{t).v(ng)'ir\t) (4.10) 

with U(t) = W{x ( R \x R ) ■ V(Sl x /). 

Therefore, the operator K( t )(n») has an iso-spectral evolution in time, and thus, its 

eigenvalues, or equivalently Tr ( V XR (Qio ) ) , are conserved in time. Using Yang-Mills integral 
equation one can write V x t (0^) as a volume or a surface ordered integral: 

V^) = P 2 e*^**^^)^ = P 3 J^ d ^ VJV -\ (4.11) 

X R 

The fact that this operator is independent of the parameterization of the volume guarantees 
that the conserved charges constructed from it are also independent of this parametrization. 
One can, however, argue about the dependence of our results with respect to the choice of 
the reference point. As we saw above, by changing x R to x R , V (*) changes under conjugation 

R 

with respect to W(x R ,x R ), and so its eigenvalues (conserved charges) are invariant under 
change of reference point. In addition, the reference point is placed at the border of Q Xr 
(the spatial infinity), and the field strength is supposed to vanish there, which implies that 
the gauge field is asymptotically flat, and therefore the Wilson line depends only on the 
points and not on the path joining them. The conserved charges are also gauge invariant, 
since under a gauge transformation we have seen that V^ct)($7^> ) — > 9RV x (t)(£l'£!)g R 1 , with g R 

being the group element, performing the gauge transformation, at x^ R . Note in addition that 
if Vw(nS) has an iso-spectral evolution so does g c V (t)($l^), with g c and element of the 

X R X R 

centre Z(G) of the gauge group G. That fact has to do with the freedom we have to choose 
the integration constants to lie in Z(G), without spoiling the gauge covariance of (4.5). 

Note that there are several integral and loop space formulations of Yang-Mills theories [8] . 
Our approach shares some of the ideas of those formulations in the sense of using ordered 
integrals of the gauge potential and of the field strength tensor. However, it differs in an 
essential way because it is based on the new eq. (4.5). In addition, it leads in a quite novel 
and direct way to gauge invariant conserved charges. 
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4.2 The conserved charges for Dyons 



Consider a theory where the gauge group G is spontaneously broken into a Higgs field in the 
adjoint representation. Consider a BPS dyon solution: 

Ei = -F 0i = sin 9 Did) B t = - -e ijk F ik = cos 9 D i( b 

2 

with 6 an arbitrary constant angle. 

At spatial infinity (r — » oo, f 2 = 1), we have 

where G(r) is an element of the Lie algebra of H, covariantly constant D^G(r) = 0[9]. Also, 
the field strength goes to zero there, which leads to a asymptotically flat gauge field 

A„ = '-d^WW- 1 

and therefore, at <S^, we get 

G(r) = WGrW- 1 

where Gr stands for the value of G(f) at xr. 

Then, the operator V(fl) discussed previously becomes 

V (t)(Q,®) = P 2 e e ^ s ™ {t) dTda \ aF ™ + ^™^^^ F ' = e [-»e(a cos6»+/3 sin8)G R ] 
x R 00 

and the conserved charges are given by the eigenvalues of Gr, which contain the magnetic 
and electric charges of the dyon solution. It is important to remark that since Gr is at spatial 
infinity, according to our observations a change of the reference point makes it change by 
conjugation with the Wilson line of the path between the new point and the old one; this 
changes nothing in the eigenvalues of Gr. 
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Appendices 



We give in the next two appendices the proofs, following the arguments of [1, 2], of the 
standard non-abelian Stokes theorem for a 1-form connection and then its generalization for 
a 2-form connection. 



A The "standard" non-abelian Stokes theorem 

Consider a Lie algebra valued 1-form C = C ll (x)dx >l , defined in a d + 1 dimensional simply 
connected space-time M. Given a path 7 parametrized by a G [0, 2ir], such that = 0) = 
xr and x fl (a = 2ir) = Xf, the Wilson line is constructed from 4 

^ J — iy 7 [a,0] =0, (A.l) 

with initial condition W[0, 0] = Wr. 

After integrated, equation (A.l) gives the path-ordered integral 



/ r dx^ \ 
W>,0] = P x exp I- J C^—da'j ■ W R . (A.2) 



There are many paths one can choose to link xr to Xf, and a natural question is whether the 
Wilson line depends upon this choice. Being M simply connected, a different path 7' with 
the same end points xr and Xf is related to 7 by continuous transformations and the Wilson 
line reads Wy^TT, 0]. 

Lets introduce a new parameter, r e [0, 2tt], labeling the path one chooses to go from xr 
to Xf. For 7 we set r = 0, and for 7', r = 2ir. All other values give intermediary paths 
between these two, that arise due to the variation process 7 — > 7 + 5-y. With that, each point 
is now characterized by two parameters, r and a: the first tells in which curve the point is, 
while the second, where in this curve. 

The variation is performed as follows. At a given point, immediately after xr we define a 
vector T M = joining 7 and 7 + 57. Basically we are defining a way to map one point of the 
curve at r = to a point of the curve at r 7^ 0. In order to answer if this holds for every other 
point, we need to parallel transport the vector along 7, i.e., along the direction of = 
the Lie derivative of in the direction of gives £ 5 T" = S u d u T»-T v d u S» = = 

and therefore if we define the variation to be orthogonal to the curve 7 at a given point, this 



4 The notation Wj[a, 0] says that the quantity W is evaluated on the path 7 from the point with a = to 
the point x M (er). 
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Figure 6: The variation of a path with fixed end points. 



will be the case for every other point. A variation of the path implies a variation of the 
holonomy, which can be calculated taking the variation of the equation (A.l): 



, ,'dW~ . s dx>* TIZ , 



Computing this explicitly 5 (omitting some symbols that we shall reintroduce at the end): 

W' 1 ^(5W) + W- 1 C, d ^5W + W- l s(c„ d ^-)w = 
da da \ da J 

(w~ 1 sw) + w~H (cj^j w = 



d 

da 



d 

da 



where after the variation was performed we multiplied the equation by W~ l from the left 
(line 1), used the chain rule in order to rewrite the first term (line 2), and with the identity 
dV ^ r7 1 = —W^^rW -1 and equation (A.l), two terms were mutually canceled (line 3). This 
last equation can be integrated from x M (cr = 0) to x^(a), a point of 7. Taking into account 
the fact that the holonomy does not change at the initial point (5Wq = 0), we get 

r f dx fl \ 

5W 1 [<t,0] = -W 1 [<t,0] J da , W- 1 [a',0)5\C l ,—jW J [a , ,0}. 



5 We introduced a new, but obvious, notation: <5W 7 [t; a, 0]. This stands for the variation of the Wilson line 
at a point x^{a) in the curve r. 
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The final step is to compute the variation appearing inside the integral above: 



da' W-H [ I W 



W-\a)C^ (x {a)) W(a)5x»(a) + / da' \ W^SC^W - — (W^C^W) 8x" 

J q l da da 

W~ l C^W5x^ + J da' W- 1 id v C„6x v — 



, dx^ dC„ 



5x> 1 + 



dx 



C ^ Cu da> 



,>V*j. W = 
SxA W 



C a f dr^ dr v 

W- 1 C li WSx' i + / da' W~ x \ d v C^x v — - d v C»—bx» 
J Q [ da da 

f a dr^ 
W- l (a)C„ {x (a)) W{a)5x»{a) - / da' W- x F. v W—bx v 

Jo da' 

where = d^Cy — dyC^ + [C^, C„] is the curvature of C M . 

From line 1 to line 2 we performed an integration by parts and used the fact that 5x^(0) = 0. 
Then, since the variation in the connection is due to a variation in the space-time point 
(x — > x + Sx), we took 5C M = d u C^5x u in line 3, and performed the derivative. After that, 
we used chain rule for the derivative of the connection in a', and equation (A.l) inside the 
commutator, which gives line 4. Then, relabeling the indices we get a curvature, showed in 
line 5. 

Finally, the variation of the Wilson line due to a variation of the path 7 to another path, 
labeled by r reads 

8W 7 [t]<t,0] = -C ti W 1 [r;a,0]5x fl (a) + 

r dx^ 

W^[r;a,0}J da' W~ l [r- a\ 0]F^W 7 [r; a\ 0]— 5x v (A.3) 

Notice that this is a general result for the variation is not specified: one can do it in both 
directions, tangent to the path, or orthogonal to it. Of course, in the tangent direction, a 
variation is just a reparametrization a — > a = f(a) , which, from equation (A.l), changes 
nothing as long as f(a) is monotonic. On the other hand, in the orthogonal direction one 
has 5x^ = T» and SW[t] = W[r + 5t] - W[t] = ^5t, which, in (A.3) (with a = 2vr, so 
that the end points are the same for both paths, and therefore 5x^(0) = 5x^(2ir) = 0) gives a 
differential equation for W: 



d r 2n dr^ dr u 

—W 1 [r;2n,0]-W y [r;27i,0} J da W; 1 ^; a,0}F^W,{T; a,0} — — 



0. 



(A.4) 
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What we conclude is that there are two possible ways to calculate the Wilson line of a given 
curve. Consider, for instance, the path 7'. The Wilson line Wy[r = 2tt; Xf, xr] can be 
obtained, first, from equation (A.l), after integration over 7'. The second possibility is to 
think of 7' as obtained from 7, using the variations we discussed above. Then the holonomy 
there is given after integrating (A. 4) from r = to r = 2tt: 



(/*2"7r p2n 
J dr J daW~, l F^W, 



2:7 r2:: r/.r'' dx" 

1 da dr 



where P 2 stands for the "surface-ordering", i.e., the ordering due to the non-abelian character 
of the fields, according to the way we vary r. 

It is clear that the quantity appearing in the r.h.s of the above equation is the flux of 
F^ v = W^F^W through the space-time surface £ C M, whose boundary is <9£ = V = 7 /-1 -7; 
lets call it §(F W , £), and rewrite the equation as 

W Y = W,-P 2 e^ FW ^. (A.5) 

Consider the case where the path discussed above is a loop T with no self intersections: 

: a E [0, 2tt] -»■ M a^(cr); z"(0) = x m (2tt). 

Then, integration of (A.l) for the whole loop gives 

W r = P 1 e$r c » dx * -W R . (A.6) 

Being T the boundary of a two dimensional submanifold S C M, we also have that this same 
Wilson line can be calculated taking a variation from the point loop P XR = Xr\/ a 6 [0, 2tt] 
(whose Wilson line is Wr) 

W r = W R -P 2 e^ FW V. (A.7) 

The fact that it is possible to calculate Wr integrating the connection C over T, or integrating 
the curvature F(C) over the area bounded by T is exactly the statement of the Stokes theorem: 

P l exp Qf^ CpdxTj -Wr = W r -P 2 exp (q>(F w , S)). (A.8) 



30 



B A generalization of the non-abelian Stokes theorem 
for a 2-form connection 



Let us consider the antisymmetric field 6 B^(x), denned in the d + 1 dimensional space-time 
M. A family of (homotopically equivalent) loops with base point xr can be used to scan the 
two-dimensional hypersurface E, starting from the infinitesimal loop around xr and taking 
variations along the direction, as explained in the previous section. 

A new quantity 7 V is introduced in analogy with the Wilson line in equation (A. 4), defined 

by 

' /U:r01 Vz[t,0] r^^M.OlVrM-of^ - (B.l) 



with the initial condition V[0,0] = Vs H , being the infinitesimal surface around xr. 

The surface-holonomy Vs [r, 0] is defined on the surface whose boundary is the loop labeled 
by r. The quantity 



/■27T 

T 27r (B, C,t)= / da Wf^r; a, 0]B^W r [r; a, 0] 
Jo 



dx^ dx v 
da dr 



plays the role of a "non-local connection" 8 , defined on each loop. The Wilson lines appearing 
inside the integral are calculated from (A.l), running from xr to x^(a) on the loop r. 

After integrating (B.l) one gets the surface-ordered integral 

V*[t, 0] = V r - P 2 exp Qf dr' T 2w (B, C, r'^j . (B.2) 

In analogy with the fact that the Wilson line is defined over the path that links two (boundary) 
points, the surface-holonomy is defined over the surface that links two (boundary) loops. Take 
these two loops to be those at r = and r = 2ir, fixed. Clearly one can use different surfaces 
to link them. Consider S and two possibilities. Each of these choices might lead to a 
different solution of (B.l). Since £' can be obtained from S after a variation in the = ^ 
direction 9 ( ( G [0, 2vr] ) we can calculate the difference 5V = V[( + 5(] - V[(] = in 
analogy with what we did for the path-holonomy, taking first the variation of the defining 
equation (B.l): 



6 The 2-form B = ~B IMV (x)dx^ 1 A dx v is not necessarily exact. 
7 For future reference we call it surface-holonomy. 

8 It was shown in [1, 2] that it is in fact a 1-form connection in the loop space. 
9 This direction is in the normal direction to the surface S. 
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Making it explicitly: 



dr 



dV 



(6VV- 1 ) - 6V ( — + T 2n (B, C, r)V- 1 ) - V6T 2n (B, C, r)V- 1 



d 

'dr 



(8VV- 1 ) -VST^BtC,^- 1 = 



We took the variation, and multiplied by V 1 from the right. After that the chain rule was 



used to get the first two terms in line 1. Then, using the identity , = —V -^V and 
equation (B.l) the term inside the parenthesis in line 2 vanishes. The next step is to expand 
the term VT 2 - k (B,C,t)V~ 1 appearing above. 



2- 



VST 27T (B, C, t)V~ 1 = V I da {5W~ l B^W + WB^SW' 1 + W6B flu W~ 1 ) -^--j-V' 1 

J o CL(T &T 

f 2n , dSx^dx" , f 2n , dx fl d5x u , 
+ V I da W- l B^W— —V + V / da W'^BJW- —V~ x 



da dr 



p2n 

V da{ [B™, W^SW] + W- l d p B^W5x p ) 
Jo 



^ da dr 

dx p dx v . 



da dr 



-v- 1 



-27T 

V I da 



o 



B w 
^ da 



dx v 
dr 



+ W^dpB^W 



dx p dx v 
da dr 



+ B 



w d 2 x v \ 
^dadr) 

dv- 1 



Sx^V- 1 



+ ^~ {VT 2lT (B, C, — ^T 27r (-B, C, 5)V~ 1 — VT 2n (B, C, 5) ^ 



-27T J 

-VI da^- 
o dr 



w-'b^w^) v-Hx» 



When the integration by parts of the second term was performed we used the fact that 8x ll {a) 



dv- 



vanishes at the boundaries. For the commutators we use SW given by (A. 3) and given by 
(A.l). In line 5 we can use the identity 
in the last term using equation (B.l). 



dr 



da 

_y-id^y-i anc i we calculate the derivative 
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Plugging the result back into the equation we started with, and integrating gives 



SVv[t,0] = (VT^foC^V- 1 ) 



^s[r,0] + 



2tt 



/ (h'Yi / daW- 1 (D p B^ + D fl B l/p + D I/ B pil )W——5xf + 



+ 



f 2w dr^ dr u f 2 * dr^ 

J da [B%, T a (F, C, 5)] - J da T a (F, C, r')] ^5x» + 



+ [T 2n (B,C,5),T 2n (B,C,r')} W- 1 Vj:[r,0]. 



The second term will be called /C, so that we can write 



5Vx[t,0] = (VT 27T (B,C,5)V- 1 ) 



Vz[t,0]+KVz[t,0]. 



(B.3) 



Taking r = 2n, the first term on the RHS vanishes and we get 



SV^[2n,0] -/CVs[2tt,0] = 0. 



(B.4) 



Then, taking the variation to be along the direction so that 5V = we get the 

differential equation for V(now we drop some of the symbols that have being used so far): 



d 



V-KV = 0, 



(B.5) 



where all 5x^ in /C is replaced by an d we write it in the nicer form: 

,dx»dx u dx p 



K 



da dr d( 



p2n [ p2n 

J drVlj daW-^DpB^ + D^p + DvB^W 



+ 
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Now let us consider a closed surface E , so that the solution of the equation (B.l) reads 



2tt 



Vx a = V R -P 2 exp (J dr T 27T (B, A, r) J . (B.6) 

The surface S D is the boundary of the three-dimensional submanifold (volume) Q, so the 
surface-holonomy Vs a can be calculated from variations in the direction starting at the 
infinitesimal surface around the reference point Xr, using (B.5): 

Vx = P 3 exp (J** d( JCj V Sr . (B.7) 

The fact that there are two ways to compute the surface-holonomy Vs is the statement of 
the Stokes theorem: 

V^ R P 2 exp (J^ dr T 2n (B, A, r)) = P 3 exp d( JC^J V^ R (B.8) 

being the l.h.s computed over the surface E which is the boundary of the volume Q, where we 
integrate the r.h.s. In order to make it more explicit we rewrite this as 

V^ R P 2 exp I J Jt da W^B^W——) = P 3 exp I J d( ICjV^ R . 
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